Slow magnetic dynamics and hysteresis loops of the bulk ferromagnet Co7(Te03)4Br6 
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Magnetic dynamics of a bulk ferromagnet, a new single crystalline compound Co7(Te03)4Br6, 
was studied by ac susceptibility and the related techniques. Very large Arrhenius activation en- 
ergy of 17.2 meV (201 K) and long attempt time (2 • 10~^ s) span the full spectrum of magnetic 
dynamics inside a convenient frequency window, offering a rare opportunity for general studies of 
magnetic dynamics. Within the experimental window the ac susceptibility data build almost ide- 
ally semicircular Cole-Cole plots. Comprehensive study of experimental dynamic hysteresis loops 
of the compound is presented and interpreted within a simple thermal-activation-assisted spin lat- 
tice relaxation model for spin reversal. Quantitative agreement between the experimental results 
and the model's prediction for dynamic coercive field is achieved by assuming the central physical 
quantity, the Debye relaxation rate, to depend on frequency, as well as on the applied field strength 
and sample temperature. Cross-over between minor- to major hysteresis loops is carefully analyzed. 
Low-frequency limitations of the model, relying on domain wall pinning effects, are experimentally 
detected and appropriately discussed. 

PACS numbers: 75.60. Ch, 75.78.Fg 



I. INTRODUCTION 

Magnetic hysteresis, historically the most intensively 
investigated nonequilibrium phenomenon, is only re- 
cently getting better understood in all necessary details 
and aspects. Traditionally, the subject has a more then a 
century long history established in research and the appli- 
cations of ferromagnetic materials in bulk- and thin-film 
forms. In the latter systems hysteresis relies on dynam- 
ics of pinned domain walls^^ a subject playing a pivotal 
role in various fundamental aspects of magnetisnP and 
its applied outcomes^J Nowadays, the subject expands 
to nanomagnetic systems^, i.e., to systems of magnetic 
nanoparticlea^, as well as to molecular- and single chain 
magnets^, not necessarily involving magnetic long range 
order at all^. Besides fundamental reasons profound un- 
derstanding of hysteresis and the related magnetic dy- 
namics is urged by the demands of rapidly growing ap- 
plications varying from, e.g., the ultra-high density mag- 
netic recording to magnetic hyperthermia in anti-cancer 
medical therap3p'. 

Let us briefly recapitulate the main physical ingredi- 
ents involved in the problem of magnetic hysteresis and 
the underlined dynamics^. A central phenomenon that 
any magnetic hysteresis relies on is a dissipative spin 
reversal. Although an atomic-scale event, spin rever- 
sal is a subject to hierarchy of interactions rendering 
the hysteresis problem mesoscopic in its character. In 
their hierarchy, interactions generate both the equilib- 
rium and non-equilibrium conditions for reversals. Equi- 
librium one follows from the appropriate free energy func- 
tional (known as micromagnetic free energ}/^) listing the 
relevant energies: exchange interaction, microcrystalline 
anisotropy, sample- geometry-dependent magnetostatic-, 
and the Zeeman energy. Minimization of micromagnetic 



functional defines then, depending on strength of the ap- 
plied field, the equilibrium single- or multi-domain state. 
The archetypal Stone- Wohlfarth hysteresis model^, e.g., 
follows from such an equilibrium considerations of a 
single-domain sample. 

Non-equilibrium aspects in spin reversal are the con- 
sequences of energy barriers separating the metastable 
spin states in neighboring local minimsP^l Modeling the 
free energy of a hysteretic system as a multi-variable 
function of its magnetic degrees of freedom one sweeps, 
in the course of hysteresis cycle, over a complicated 
energy landscape. The barriers themselves rely on a 
number of possible sources, like a domain wall pinning 
on impurities/imperfections, spatial variations in local 
anisotropics, etc. Energy barriers are the key elements 
of magnetic dynamics- on basis of the Arrhenius thermal 
activatioii^*^ they introduce, e.g., the relaxation time ef- 
fects . 

Finally, spin reversal is a cooperative phenomenon- 
reversing spins belong collectively to the mesoscopic or 
macroscopic objects, domain walls, involving large num- 
ber of spins. The details are material-specific, via the 
peculiarities of the free energy landscape. 

Out of different features characterizing hysteresis the 
most prominent one is coercive field- a characteristic re- 
verse field inducing massive nucleation/propagation of 
domain walls. The ambition of any hysteresis model 
is to provide a realistic prediction for the the coercive 
field. Usually, the models report their results formulated 
as the corrections to the elementary Stoner- Wohlfarth 
result by taking into account dynamical features and 
effects^. In this article we present experimental study 
and appropriate modelings approaching the hysteresis 
problem the other way around- a central investigated 
topic is dynamical hysteresis while the hypothetical dc 



hysteresis is considered as a dynamical one taken in 
the low- frequency limit. Opened freque ncy window is 
very important in hysteresis studies^^^^ from the num- 
ber of reasons. In coercivity studies on real samples the 
frequency- sensitive methods help, for example, to dis- 
tinguish between the contributions originating from the 
Arrhenius (or Neel-Brown) thermal activation or from 
the static micromagnetism^^. There have been a lot of 
efforts invested in the attempts to reconcile the predic- 
tions of these two approaches and to identify the con- 
ditions (like temperature, frequency range) favoring one 
over the othei^^ - this work keeps on the similar track. 

In this study we focus on the low- frequency, domain 
wall-related, magne tic dy namics of a new oxo-halide sys- 
tem Co7(Te03)4Br6^^ ^^. Owing to rare combination 
of parameters determining its magneto-dynamics, con- 
veniently matching the dynamical hysteresis' experimen- 
tal window, we show that the experimental hysteresis 
can be quantitatively reproduced within a simple Debye- 
relaxation-based model for thermally activated spin re- 
versals. For this model to be successful in the interpreta- 
tion of spin dynamics, as represented by the correspond- 
ing dynamic hysteresis loops, one has to allow the relax- 
ation time to depend on frequency and field, as well as 
on temperature. 

The article is organized as follows. The involved ex- 
perimental techniques, sample information, as well as the 
main motivations, are introduced in Section|TI[ In Section 
III we show that domain wall dynamics just below the fer- 
romagnetic transition (of the first order type) obeys the 
Debye relaxational dynamics in a wide frequency range 
of the presented Cole-Cole plots. We also discuss the rea- 
sons for the observed exceptionally long relaxation times 
and elaborate why the conditions are favorable for in- 
depth studies of magnetic dynamics. Experimental in- 
sight into various stages of the spin-reversal-based do- 
main wall dynamics is presented in that section. Section 
JV|presents a simple model of thermal-activation-assisted 
Debye relaxations comparing afterwards the model pre- 
dictions with the detailed experimental study of dynamic 
hysteresis of Co7(Te03)4Br6. In Section[v|we discuss the 
low-frequency limitations of the model, drawing out ap- 
propriate conclusions. 



II. SAMPLE AND EXPERIMENTAL DATA 

The measurements were performed on single crys- 
tals of a ne w, c obalt-based oxi-halide compound 
Co7(Te03)4BrJ^2Eil. Magnetism of this system is char- 
acterized by strong magnetic anisotropy extending far 
into the paramagnetic temperature range. The mea- 
sured bulk anisotropy relies on the strong single ion 
anisotropy energy of variously coordinated Co^+ ions. 
At T/v=34 K an incommensurate antiferromagnetic order 
sets in. Its propagation vector is temperature dependent 
and, at Tc=27.2 K, an abrupt spin reorientation transi- 
tion introduces a ferromagnetic component along 6-axis. 



Standard phenomenology of ferromagnetism character- 
izes measurements along this axis onl};^^ -along other 
axes standard phenomenology of antiferromagnetic or- 
der, like spin-fiip (or fiop) transition, persists below Tc. 

Ferromagnetic component of magnetic order of 
Co7(Te03)4Br6 sets in accompanied by extremely sharp 
and sizable peak^^ ^^ in ac susceptibility. Surprisingly, 
the peak in imaginary susceptibility {x") was found, de- 
pending on the strength and frequency of the applied 
field, significantly bigger than the peak in real suscep- 
tibility ix')- The latter result is at variance with the 
usual finding x" ^ x' dominating by far magnetic tran- 
sitions studies of various magnetic systems in general. 
In an attempt to figure out why Co7(Te03)4Br6 devi- 
ates so much from the standard behavior we realized that 
it represents in fact a rather unique system: its intrin- 
sic magnetic parameters discussed below enable the en- 
tire spectrum of magnetic dynamics to unfold within the 
most convenient /sensitive experimental window, i.e., in 
the frequency range 0.05 Hz- 1 kHz. It thus enable a 
direct experimental access to various stages of magnetic 
dynamics by ac susceptibility and the related techniques. 
A surprising observation mentioned above, x" ^ X^ ^ i"^- 
mark of big dissipation exploding in the transition, can 
be then very simply interpreted. 

The ac susceptibility/dynamic hysteresis studies were 
performed by the use of a non-SQUID CryoBIND ac sus- 
ceptibility system, revealing high (2 nanoemu) sensitiv- 
ity, within the frequency range 50 mHz-1 kHz. The sam- 
ples used were oriented single crystals cut into the rod- 
like geometry (typical size 0.1x1x4 mm^), thus minimiz- 
ing the value of the demagnetizing factor. The longest 
sample axis corresponds to the ferromagnetic 6-axis and 
all measurements were performed with dc or ac field ap- 
plied parallel to this axis. 



III. RELAXATIONAL DYNAMICS OF DOMAIN 
WALLS 

In order to clear out nature of the ferromagnetic tran- 
sition at Tc we first show in Fig. [l] a set of dynamic 
hysteresis loops of Co7(Te03)4Br6 (discussed at length 
in the next Sections). In the transition region these 
loops were taken within a narrow, 0.14 K-wide, tempera- 
ture interval. Ferromagnetic order parameter -saturation 
magnetization- builds up its low-temperature (4.2 K) 
value practically within 0.1 K, illustrating the first-order 
character of the ferromagnetic transition^^. This means 
that the magnetocrystalline anisotropy, which inherits 
its temperature dependence from the temperature de- 
pendent saturation magnetization^^, is fully developed 
already at 27 K and barely changes by lowering temper- 
ature below 27 K. Simultaneously with anisotropy and 
bulk magnetization it is clear from the first principles^ 
that magnetic domains and the domain walls set in and 
self-organize just below T^, in order to keep the magnetic 
free energy minimal ; The domain pattern then responds 
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FIG. 1. Dynamic hysteresis loops of Co7(Te03)4Br6 taken 
at several temperatures within the ferromagnetic transition. 
Numbers adjacent to the loops designate the temperatures. 
A small hystersis assymetry -a shift to the left- is commented 
in Section IIVI 
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FIG. 2. (Color online) Frequency dependence of x (open 
symbols) and x' (solid symbols) of Co7(Te03)4Br6 just be- 
low Tc — 27.2 K (main panel). Solid lines are fits to the stan- 
dard expressions for the relaxing complex susceptibilit3/^(Eq. 
(3)) with the addition of a prefactor in the formula for x'' 
to account for its asymmetry. Dynamic hysteresis and dissi- 
pation scans were measured in the tail area characterized by 
x' > X • Inset shows the same data from the main panel but 
plotted in the x^ X diagram- the data are now collapsed on 
the single Cole-Cole semicircle. 



to any changes in the parameter values by dissipative do- 
main wall rearrangements. 



A. Cole-Cole plots 

Frequency dependence of ac susceptibility, x{.^\ Pro- 
vides the most elementary insight into magnetic dynam- 



ics of ordered systems . Fig. [2] shows the frequency de- 
pendence of the susceptibility components x!^^\x!,\^) 
of Co7(Te03)4Br6 crystal taken at several tempera- 
tures below Tc. This plot provides a detailed insight 
into magnetic dynamics of the compound. The type 
of response shown in Fig. |2] tells that the responsi- 
ble magnetic subsystem, in its response to the applied 
field, interacts strongly with the lattice: the pattern 
of x'{^iT)^x"{^iT) curves is reminiscent of a sim- 
ple Debye-type spin lattice relaxation^. Indeed, the 
x'{oJ^T),^{uj, T) data, if presented in the form of Cole- 
Cole plot^^, i.e., as the functional form x" {x') (Inset), 
all collapse on a somewhat distorted Debye semicircle 
(common for all temperatures within the range^. Ideal 
semicircle, a fingerprint^^ of single relaxation time r, 
intersects the horizontal x'-ci^xis in the values of adia- 
batic xs (=0? A' ^ co) and isothermal xt (^ ^ 0) 
susceptibilit}'^^. 

A closer examination of x' {^) •> x" i'^) data, as well as 
the fact of saturation magnetization relaxing logarith- 
mically in time (data not shown), is consistent with a 
distribution of relaxation times. However, in order to 
grasp only the basic mechanisms we hereby focus just 
the characteristic relaxation time r, as determined by the 
frequency maximizing x" • ^ niay be considered as the av- 
erage time the magnetic subsystem, taken out of equilib- 
rium by the application of, e.g., external magnetic field, 
requires to thermalize again. The temperature depen- 
dence of r (Fig. [3| is obviously subject to the Arrhenius 
thermal activation. The involved activation energy Ea 
represents the effective energy barrier for the thermally 
activated spin reversal. The reversal applies collectively 
to spins belonging to certain volume Vq of the domain 
wall, an elementary reversal entity. In the context of 
domain wall dyna mics Vq is usually termed Barkhausen 
or switching volum^^^ while in superparamagnetism- 
related literature the related term activation volume is 
more common. (Although there are differences in precise 
physical meanings of these termP^ these differences are 
not of big importance for the main focus of this article) . 

The elementary switching volume Vq for 
Co7(Te03)4Br6 can be readily estimated from 
Ea = ^Vb, where K represents the magnetocrys- 
talline anisotropy energy. Estimating crudely K from 
the Stoner-Wohlfarth approximation. He = 2K/M, 
and the measured DC susceptibility values^^ for He 
and M one gets Vq = 5.5 • 10~^^ m^. The edge length 
of the corresponding cube equals 38 nm, comprising 
approximately 1500 Co-ions. Activation volumes of very 
similar sizes have been frequently reported for, e.g., 
ultrathin Co-based magnetic films^. 

The plots shown in Fig. |2] witness that magnetic dy- 
namics of Co7(Te03)4Br6 is exceptionally well exposed 
to experimental studies. A conventional ac susceptibility 
frequency window, 50 mHz - 1 kHz, is usually too nar- 
row to span the dynamics in its full spectrum. This is 
reflected by a narrow arc spanned in the x^^ix') plots of 
many investigated magnetic systems, instead of the ideal 
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FIG. 3. Temperature dependence of the mean relaxation time 
(main panel) obeys the Arrhenius-like thermal activation with 
approximately temperature independent activation energy Ea 
(inset). Ea is unusually big due to strong magnetocrystalline 
anisotropy K resulting from the strong single- ion anisotropy 
energy. 




FIG. 4. (Color online) Plot of the Arrhenius thermal acti- 
vation T{T,Ea,ro) = Toexp{Ea/kBT) at fixed temperature 
T=26 K, in its dependence on the attempt time tq and the 
activation energy Ea. Thick black curves defines approxi- 
mately the experimental window involved with the present ac 
susceptibility study. Intrinsic parameters of Co7(Te03)4Br6 
(values of Ea and tq) at T=26 K impose magnetic dynamics 
to take place in the position of the black dot, i.e., by coinci- 
dence, within the experimental window of readily observable 
area on the In t(to, Ea) surface. By further cooling magnetic 
dynamics slows down -the ln(T) surface rapidly elevates- and 
the system's dynamics status -black dot- rapidly escapes from 
the experimental window area. 



Cole-Cole semicircle^. In Co7(Te03)4Br6, however, the 
complete spectrum of the low- to high-frequency modes 
gets activated within a conventional frequency window. 
This situation represents obviously a consequence of un- 
usually slow relaxations characterizing magnetic dynam- 
ics of Co7(Te03)4Br6, established already 1-2 K below 
Tc. Its slow relaxations are represented by long relaxation 
times r{T) at the order of second, Fig.lsj determining the 
Arrhenius activation quantitatively. 

Physically, slow relaxations stem from the large values 
of both Ea {Ea/kB = 201 K) and tq (= 2-10-'^ s). Let us 
illustrate schematically the experimental significance of 
these figures for studies of dynamics. Fig.[4]shows the Ar- 
rhenius surface r{ro^Ea)^ as it looks like on logarithmic 
scale for T=26 K, and the position of Co7(Te03)4Br6, 
as determined by its intrinsic values of tq , £^a • If these 
values were much smaller, as they usually are for large 
number of ferromagnets studied thus far, the position of 
the system at the 26 K-r(ro, Ea) surface would be deep in 
the blue/dark zone, i.e., well away from the experimen- 
tally achievable frequency window. With actual values of 
To^Ea for Co7(Te03)4Br6 its dynamic status stays inside 
the experimental window for temperatures in a pretty 
broad range, 22 K - 27 K, enabling detailed studies of 
domain wall dynamics. 

We relate a sizable Ea primarily to the fact that 
ferromagnetism arises in this compound out of already 
macroscopically anisotropic, AF ordered, backbone^^. In 
this pre-established order a strong single-ion anisotropy 
plays important role, as well as in the microcrystalline 
anisotropy energy Ea setting-in abruptly below Tc. As 
for the attempt time tq its experimental value is orders 
of magnitude bigger than the microscopic spin-flop at- 
tempt time at the order of 10~^ s. Big tq is naturally 
attributed to collective dynamics of dom ain w all spins, 
known to take place on a slower time scale^^I^. 



B. Dissipation scans 

With the relaxation time conveniently positioned 
within the experimental ac susceptibility window there 
is an additional freedom in choosing a particular ap- 
plied field's working frequency in the relaxation diagram. 
Fig. [2] For examples, one can choose frequency on the 
right-hand side of the peak in x"{f) (consistent with 
x" ^ X^ i-6., with high dissipation) and do the tempera- 
ture or the dc-field scans. With such choice of frequency 
the temperature scan reveals explosive peak in x" {"E) at 
Tc mentioned in the Introduction Section, which can now 
be reproduced by quantitative precision within a simple 
model based on the expressions for x' iX" i^ the Debye 
approximation^ (data not shown). More interesting are 
the 'dissipation scans' taken with the same choice of ap- 
plied frequency. Dissipation scans hereby refer to record- 
ing of x^', a measure of dissipative magnetic dynamics, by 
the use of a small ac field amplitude superimposed over 
the dc field slowly ramping up and down. Fig. [5] DC field 
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FIG. 5. Dissipation scans (see text), shown at different tem- 
peratures, measure the rate of spin reversals as the ramping 
dc field progressively ejects, and then nucleates back, the do- 
main walls. Field sequence is indicated by horizontal arrows 
in (a). The points ascribed to domain walls' nucleation (N) 
and propagation (P) are indicated by the inclined arrows. A 
smooth background dissipation is attributed to incomplete 
spin reversal. The virgin curve dissipations are shown in (d). 



domain wall dynamics, otherwise taking place in stan- 
dard dc hysteresis cycling. 

In the scan-starting condition of a zero-field-cooled 
sample at Hdc = 0, necessarily composed of equal num- 
ber of -\-M and —M domains, the measured dissipation 
was already at a high initial level. Fig. ^). To achieve 
the highest possible rate of spin reversals application of 
a small temperature-dependent dc field is obviously re- 
quired (Fig. [sji)). We speculate that a maximum in dis- 
sipation is achieved after the field supplies enough mag- 
netic energy to overcome a weak domain wall pinning 
energy. In the virgin curve a rapid dissipation decay 
indicates ejection of the domain walls out of the sam- 
ple rendering it single-domain and revealing just residual 
dissipation afterwards. In the field back-swing reversed 
domains start nucleating again and further dissipation 
evolves through propagation of the domain walls as well. 
Both processes are integrated into the 'dissipation dome' 
structure. Fig. [5] The temperature-dependent positions 
of the maxima can naturally be identified with the coer- 
cive field Hc{T). The Hc{T) values, being small immedi- 
ately below Tc but well-defined, are in good quantitative 
agreement with the dc-SQUID Hc{T) datsP. 

The structure of the dissipation dome tells not only on 
the activation of the dissipative nucleation/propagation 
modes but also on the importance of thermal activa- 
tion. After being expelled from the sample the domain 
re-nucleation relies on overcoming the anisotropy-related 
energy barriers for spin reversal of the activation vol- 
ume. This obviously takes place with the most frequent 
rate at He revealing, naturally, the field strengths be- 
ing negative in sign (in the field first back-swing). But 
one also notes, for loops taken at T=26.5 K, 24 K, that 
the very onset of nucleation takes place already for posi- 
tive field values. Domain nucleation in positive field can 
only rely on thermally-induced reversals of the activation 
volumes. The latter finding provides an independent ar- 
gument favoring the importance of thermal activation in 
studies of magnetic dynamics in Co7(Te03)4Br6. Fig. p] 
also shows that, unlike the domain walls ejection from the 
virgin state, the il^c-build up is not a smooth process. Ob- 
viously, at lower temperatures it progresses via erratic, 
perhaps avalanche-based, mechanism which is not well 
understood as yet (and is not elaborated in this work 
any further). 



IV. DYNAMIC HYSTERESIS LOOPS OF 
Co7(Te03)4Br6 



sweeping is accompanied by continuously changing pat- 
tern of magnetic domains through processes of domain 
nucleation and/or domain wall propagation. On atomic 
scale these processes are underlined by dissipative spin 
reversals. Thus the dissipative scans tell directly about 
the average rate of spin reversals in successive stages of 



Now we present a more detailed study of the evolu- 
tion of magnetic dynamics in Co7(Te03)4Br6. On e way 
to do it, as shown in the studies on ultrathin films^^^, 
would consist of collecting a set of the Cole-Cole plots 
for different Hac- In owe work, however, the magnetic 
dynamics was studied more directly, by measuring the 
induction-type dynamic hysteresis responding to the si- 
nusoidal applied field Hac = ^o sin ujt. The M — H loops 
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FIG. 6. An example of dynamic hysteresis measurement on 
Co7(Te03)4Br6 in the time domain. Magnetic field pulse 
Hacif) and the induced voltage Vindit) (circles) are measured 
and recorded directly while magnetization {M{t) (triangles) 
is obtained by numerical integration of Vind{t)- The phase 
shift between Hac{t) and M(t) brings a direct information on 
the involved magnetic dynamics (see, text). All quantities are 
plotted scaled by their maximal values. 



were determined at different temperatures and studied 
as functions of the field amplitude Hq and frequency uj. 

Magnetization M{t) of the sample was determined by 
fast digital acquisition of voltage Vind{t)^ induced in the 
highly balanced secondary coil of our ac susceptometer by 
the application of sinusoidal current through the primary. 
By virtue of Faraday law, Vind{t) ^ d(j)/dt ~ dM/dt. 
Fig. [6] shows the relevant signals in the time domain. 

The period of the applied pulse was varied within the 
range 1 ms-1 s enabling studies of magnetic dynamics 
on the 1 Hz-1 kHz frequency scale. In order to improve 
statistics the identical subsequent strings of Vindif) data 
were appropriately averaged. Depending on the field fre- 
quency averaging takes 1-100 s, enabling continuous mea- 
surements of appropriately averaged Vind{t) data in the 
course of a slow temperature drift. The field signal Hacif) 
was taken by recording voltage on the standard resistor 
connected in series with the primary coil. In numerical 
elaboration of the data, taken afterwards, magnetization 
M(= J Vind{t)dt) was obtained by numerical integration 
and subtraction of the background. M{t) is then plotted 
versus Hac{t) to get the hysteresis loops in their temper- 
ature evolution. 

In Fig.[7|a set of representative dynamic hysteresis data 
on Co7(Te03)4Br6 is shown. At fixed chosen frequencies 
and measurement temperatures a sequence of loops were 
taken in increasing magnetic field amplitudes, defining a 
broad range of conditions for dissipative spin reversals. 

Referring to the loop dependence on Hq there is ob- 
viously a crossover between the two generic loop types, 
i.e., between the minor loops for small values of Hq and, 
in increasing Hq^ the major loops revealing characteris- 




-150 -100 -50 50 100 150 -150 -100 -50 50 100 150 



0.04 
0.02 


. 23 K 
29 Hz 




c) - 








0.00 


^^^ 


-0.02 
-0.04 







24 K<T<27.5 K 
29 HZ/.' ^^^''- 




-150 -100 -50 50 100 150 -60 -40 -20 20 40 60 
Hdc (Oe) H,, (Oe) 

FIG. 7. (Color online) Experimental, induction- type dynamic 
hysteresis on a 2.9 mg single crystal taken with sinusoidal 
magnetic field pulses varying in amplitudes (5-100 Oe RMS). 
Calibration to absolute units was performed by the use of dc 
susceptibility results. The pulse widths were t =34 ms and 
t =0.59 s, shown in panels (a),(c),(d) and (b), respectively. 
At higher frequencies (29 Hz) the loops -generic ellipses- are 
characteristically horizontally aligned (panel c)) reflecting ful- 
fillment of the condition x' ^ X (see, text). At lower tem- 
peratures the domain wall propagation does not reach the 
sample boundaries within the available field amplitudes. All 
aspects of magnetic dynamics are revealed in the set of 29 Hz 
hysteresis loops at different temperatures, shown in (d), in- 
cluding the collapse of order parameter for T passing through 
Tc. The loops are also characterized by small asymmetry 
commented in the text. 



tic 'whiskers'. As it is generally well understood^ minor 
loops correspond to domain wall dynamics taking place 
in the limit of infinite sample size while the major ones 
correspond to saturation to the sample-size-limited total 
magnetization. In the whisker-like hysteresis regions do- 
main walls are entirely ejected out of the sample thus any 
dynamics dimes out in the whisker region. During the cy- 
cle dynamics reappears again within the reversed mag- 
netic field swing. It is then underlined by back-nucleation 
and propagation of domain walls. The activation ener- 
gies involved with minor and major loops are not identi- 
cal; while the domain wall propagation is mostly involved 
with the minor loops the major loops relies, besides prop- 
agation, on the domain wall nucleatioiP^'^. Thus the dy- 
namical features of minor and major loops are expected 
to be at least quantitatively different. The details of these 
differences, showing up in a continuous transformation of 
minor- into major loops in increasing Hq^ are presented 
and analyzed in this study. 

Hereby, we pay particular attention on the phe- 
nomenology of minor loops, attracting recently a pro- 
nounced renewed interestP^^. As illustrated in Fig.[7|(c) 
and elaborated below the minor loop is a generic ellipsis, 
the parameters of which depend sensitively on the applied 



field magnitude Hq, its frequency /, and the sample tem- 
perature T. To the best of our knowledge the detailed 
elaboration of the latter aspects of magnetic dynamics 
has not been presented so far for any ferromagnetic sys- 
tem. For this purpose we first introduce a simple model 
for minor loops. Then, in the next section, we show the 
model to be in sound qualitative and quantitative agree- 
ment with the presented experimental investigations of 
magnetic dynamics of Co7(Te03)4Br6. 



A. Model for minor hysteresis loops 

For the sake of interpretation of the experimental loops 
we first invoke the simplest phenomenological concept of 
hysteretic behavior in general: macroscopic 'response' of 
a hysteretic system (i.e., magnetization in our case) lags 
behind the 'external probe' (i.e., instantaneous value of 
the applied magnetic field) for the phase shift 0. Thus 
the magnetization M(t), lagging behind the field H{t), 
reads 

H{t) = Hq sin ujt^ 
M{t) = Mo sm{ujt - 4>) = Mo cos0 sin ujt - 

— Mo sincj) cos ujt. (1) 

If one then defines the in-phase and the out-of-phase 
susceptibility components as x' — i^o/Ho) cos (p and 
y^ = (Mo/i^o) sin^, respectively, the magnetization be- 
comes 



M{t) = x'Ho sin ujt - x" Ho cos ut. 



(2) 



The phase shift is now determined as the ratio of the two 
susceptibility components, tan^ = x" /x' - 

At this stage it is instructive to visualize several char- 
acteristic M{H) plots for several chosen values of (/). In 
Inset to Fig. [s] we show the corresponding M(t)^H{t) 
values (for t treated as an implicit variable in Eq.Q), 
revealing the hysteretic M — H loops. The loops are 
represented by generic ellipses tilted with respect to the 
i!f-axis for an angle determined by the phase factor cj). 

The latter simple model for hysteresis loops is very 
general as no assumption for the involved type of magne- 
tization (irrespective if it is of paramagnetic, ferromag- 
netic, superparamagnetic,..., origin), nor for the involved 
phase-lagging mechanism, has been introduced. It how- 
ever strictly applies only to the case of linear magnetic 
response, i.e., to Mq being independent on Ho^ a condi- 
tion ascribed to minor loop regime. We focus our elab- 
oration now to the specific dynamics of domain walls 
in Co7(Te03)4Br6. As shown in Section 



III 



the latter 

dynamics fits nicely the model of Arrhenius-Brown-Neel 
thermal activation. Within the latter model the dissipa- 
tive spin reversal is thermally activated thus it is natu- 
rally to assume that thermal activation plays important 
role in hysteresis formation as well; hereby, we think of 
dc-hysteresis as the dc-limits of the frequency-dependent 
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FIG. 8. Coercive field factor /i(cj,t) = ujt/\/1 + cj^r^ for 
three representative values of r, effective in thermal activation 
of Co7(Te03)4Br6 at approximately 24 K, 26.5 K and 27 K, 
have been selected for the plots. Inset: Magnetization vs. 
applied field from Eq.(IT]), normalized to respective amplitude 
values Hq , Mo , for several chosen values of the phase shift 
parameter cj) (numbers adjacent to each loop). 



dynamic hysteresis. The general model of dynamic hys- 
teresis becomes specific once the phase parameter (j) gets 
defined. We therefore introduce the phase relationship 
inherent to the domain-wall-related relaxational dynam- 
ics. In its simplest Debye form (single relaxation time r) 
and within the Casimir-du Pre approximation this dy- 
namics is represented by the well-known non-resonant 
magnetic susceptibility forni^ : 



X =X5 + 



XT -Xs 



X 



{Xt - Xs)^r 



1 



(3) 



In this simple linear response theory approach all physics 
of the system is condensed in the phenomenological con- 
stant r, which, in turn, can depend on a number of mi- 
croscopic and macroscopic parameters. Remaining con- 
sistently at the phenomenological level one can put now 
the phase shift ^, a quantity decisive for the shape of gen- 
eral hysteresis (Inset to Fig. [8|, on physical grounds by 
relating it to the relaxation time r: Taking for adiabatic 
susceptibility Xs = ^ (in accordance with our experimen- 
tal data) one gets^^ from Eq.3 that the Debye- Casimir- 
du Pres relaxation implies tan^ = x" /x' = ^'^i thus 
(j) = arctancjT. 

One naturally asks to what extent the experimental 
details of dynamical hysteresis of Co7(Te03)4Br6 shown 
in Fig. [7| can be described by such a model. As a de- 
manding test one may look for coercive field and its de- 
pendence on the strength and the frequency of the ap- 
plied magnetic fields. In this study, the term 'coercive 
field' is used strictly in the context of dynamical hystere- 
sis loops- the corresponding values differ a lot from the 
'technical' coercive fields, usually extracted from the dc 



magnetic hysteresis on samples reaching magnetic satu- 
ration. Similarly as for the dc hysteresis we define the 
(dynamic) coercive field He simply as M{Hc) =0. As 
one reads from Eq.Q magnetization extinguishes when- 
ever ut = (j), i.e., disregarding the periodicity, at instant 
tc = 0/a;. The corresponding field value (i.e., coercive 
field) is: 



H{tc) = Hc{uJ^ t) = Ho sin utc = Hq sini 

(jJT 

= Hq sin arctan ujt = Hq sin arcsin 



UJT 



Ho 



\/1+w2t2 



•(4) 



Eq. (4) represents the operative result of this simple 
model which can be tested now on the experimental dy- 
namic hysteresis curves. 

Depending on the frequency of the applied field there 
are obviously two limits. In the limit of a long measure- 
ment period with respect to the relaxation time (hys- 
teresis at low frequencies, cjr <C 1) Eq. (4) implies 
He ~ Hquot. Hysteresis axis is then approximately 
aligned along the line M{H) = MqH in the M - H 
plot, cf. Inset to Fig. |8] In the opposite limit of short 
measurement period (high frequencies, cor :^ 1) Eq. (4) 
implies He ~ Hq and the phase lagging approaches its 
maximum value (p ~ 7r/2; hysteresis axis is then hori- 
zontally aligned, cf. Inset to Fig. [8J The main panel 
of Fig. [8] shows the frequ ency dependence of the factor 
h{uj^ t) = cjr/Vl + uj'^r'^j which rules, at fixed r, the fre- 
quency dependence of He in this model (Eq.Q). 

Invoking now our experimental dynamic hysteresis 
(Fig. [t]) one easily notes that there is at least qualita- 
tive accordance: A change in hysteresis alignment, from 
a steep almost vertical one (panel (b)) to the horizon- 
tally aligned one (panel (c)), can naturally be interpreted 
as transformation in measuring conditions introducing a 
cross-over between the low- to high-frequency limits in 
the presented model. 

Summarizing the physical ingredients built-in in this 
dynamic hysteresis model we note that these ingredients 
are identical to those underlying the standard spin-lattice 
relaxation modeP for x^ x'^ usually presented in the 
form of the Cole-Cole plots (see, section [Til A ) . Thus dy- 
namic hysteresis studies, interpreted within this model, 
could represent a method alternate and/or complemen- 
tary to the usual spin-lattice x^ X^'-^'oute. For example, 
Eq. (4) offers, in its low- frequency limit, the relaxation 
rate r to be extracted from the measurements of He. 
Especially in measurements looking for temperature de- 
pendence of r dynamic hysteresis route might turn out 
to be advantageous. 

Referring again to all our experimental dynamical hys- 
teresis (Figs. 1,7,10) one notes a small but systematic 
shift of the hysteresis body (at the amount of a few- 
Oersted) to the left. Thus, if one defines H~ and H^ as 
the coercive fields for the field-down and the field-up hys- 
teresis branches, respectively, one notes | H~ \>\ H^ |. 



One reason for the loop asymmetry could be a b-axis 
component of the Earth field, not compensated in the 
presented experiments. Hysteresis shift of typically 3- 
5 Oe would however be difficult to interpret on basis 
of the interference with the Earth field (approximately 
0.3 Oe). Another possible reason for the loop asym- 
metry might rely on the fact that a subject b-axis fer- 
romagnetism in Co7(Te03)4Br6 represents just a com- 
ponent of the global antiferromagnetic order. We thus 
speculate that a scenario similar to the one valid in the 
exchange-biased structures, revealing characteristic hys- 
teresis shifts, might be effective. More work would how- 
ever be necessary to clarify the observed loop asymmetry. 



B. Experimental dynamic hysteresis loops of 

Co7(Te03)4Br6: Dependence on frequency and 

amplitude of measuring fields 

In view of the model presented above, particularly its 
Eq. (4), we present now a more quantitative study of the 
dynamic hysteresis results. Eq. (4) is obviously symmet- 
ric with respect to r and uj - this model's coercive field is 
a function of the product ujt product only. For the fixed 
UJ one can first check the coercive field dependence on re- 
laxation time r. In Co7(Te03)4Br6 r is strongly temper- 
ature dependent and, below Te^ obeys the Arrhenius acti- 
vation law r{T) = roexp{Ea/kBT), with tq = 2 x 10~^s 
and Ea/kB =201 K. Within a temperature window of 
several Kelvins the relaxation time r grows by cooling 
(Fig. |3| from zero (at Te) to several seconds. Thus for 
the applied field in the typical frequency range 1 Hz - 1 
kHz the high frequency limit u ^ 1/r will in any case be 
achieved soon below Te. In view of Eq. (4) one expects 
steep increase of He below Te followed by a saturation of 
He at Hq^ the amplitude of the applied sinusoidal field. 
This is exactly what has been observed for one set of 
measuring parameters (/, i^o). Fig. [9J Saturation of He- 
values at the level of the field amplitude Hq is accom- 
panied by leveling of the hysteresis-ellipsis axis with the 
H-axis (cf. Fig. JTJi) ), just as in the elementary model of 
dynamic hysteresis in the limit (j) -^ 7r/2 (Inset to Fig.lsl). 



1. Frequency dependence of He 

For the sake of a more detailed verification of the dy- 
namical susceptibility model a large number of the mea- 
suring sequences were taken at several stabilized temper- 
atures below Te. A measurement sequence was typically 
taken on a zero-field-cooled sample at fixed applied field 
frequency, performing acquisition of hysteresis loops us- 
ing gradually increasing field strengths. Before taking 
another sequence, at different frequency, the sample was 
routinely heated above Te and cooled back to the temper- 
ature of the previous sequence. Fig flO] shows the two rep- 
resentative sets measured at two frequencies in a broad 
range of applied fields Hq. One immediately notes that in 
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FIG. 9. Temperature dependence of He as determined from 
the dynamical hysteresis curves collected during a slow tem- 
perature drift (like those shown in Fig. u\ panel d)). Some 
sample dependence has been observed in the He (T) measured 
on different samples. Solid line is guide to the eye. 



measurements taken at 10 Hz the saturated, whisker-like, 
hysteresis were clearly achieved while at 60 Hz only the 
distorted minor loops were present. The recorded sets 
of dynamical hysteresis loops enable a quantitative data 
analysis in terms of frequency and field strength treated 
either as a variable or as a parameter. Here we present 
the coercive field dependence first on frequency and then 
on the field strength/magnitude, exploring the range of 
applicability of the model's predictions (as summarized 
by Eq. (4)). 

FigjTT] shows experimental dependence He on fre- 
quency, to be compared with the model's predictions, 
illustrated in FigjHj By increasing frequency up to 60 
Hz and above, the expected result. He = Hq^ achieved 
at all fields in the range, is clearly consistent with the 
high-frequency limit of the model. However, agreement 
within the full frequency range (solid lines in Fig 11 ) was 
achieved only after the relaxation rate r was assumed to 
represent not only a function of temperature but also a 
function of frequency and magnitude of the applied field, 
thus r = r(T, /, JT^o)- The physical background behind 
this assumption is discussed in the next section. The fit- 
ting procedure was performed in two steps. First, the 
best possible agreement between the experimental Hc{f) 
points and Eq. (4) was found assuming r to depend only 
on Hq in the form of a simple second-degree polynomial 
in Hq (thin solid lines in Fig 11). While for frequencies 
above 15-20 Hz the thin lines, representing the Hdf) val- 
ues calculated from Eq. (4), are in sound agreement with 
the experimental Hc{f) points there are large deviations 
in the low-frequency range. Frequency dependence of r 
was therefore introduced as well, by assuming a heuristic 
expression for r, r{f,Ho) = rh{Ho)rf{f). With r/(/) 
and ThiHo) both in the form of second-degree polyno- 
mial there is obviously a full agreement between the ex- 
perimental points and the prediction of the model for 
the Hc{f) dependence extending in the whole frequency 
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FIG. 10. Dynamic hysteresis on m=2.4 mg Co7(Te03)4Br6 
sample at 26.5 K in sequence of applied field amplitudes 
ifo (10-200 Oe, 10 Oe steps), at the two measuring frequen- 
cies, 10 Hz, and 60 Hz. To save space the sets of hysteresis 
loops taken at number of other frequencies are not shown. 
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FIG. 11. (Color online) Dependence of He on frequency / 
at fixed Hq. The experimental He points (circles) are de- 
rived from measurements partially shown in Fig |10| (for two 
frequencies only). Thin solid lines represent the predictions 
of the model, Eq|4] treating r as a fit variable depending on 
Ho only. Thick solid lines represent fits to Eq[4] treating r 
as a function of Hq and f in the form of a simple product, 
r(/, ifo) = Tf(f)rh(Ho). Domain wall pinning takes impor- 
tant role in the triangle area. 
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FIG. 12. Plots of the single- variable functions t/(/) and 
rh{Ho) used in heuristic expression for r, t(/, ifo) = 
Tf{f)Th{Ho). These functions are used to interpret the ex- 
perimental Hc{f) and Hc{Ho) points (FigjTT]) and FigjlS] re- 
spectively) 



range. The analytical forms for rh{Ho) and rf{f)^ found 
to give the best fit, are shown in FigjT2] We note that 
the rh{Ho) depends to a large extent on the applied Hq 
range- if a wider range were applied the fit parameters of 
Th would be very different. 



2. Dependence of He on the field magnitude: minor- 
major- hysteresis loop crossover 



to 



Focusing primarily the minor loop mechanisms we were 
focusing, up to now, the low-field regime (far from mag- 
netic saturation). Now we want to extend our results to 
higher fields. Instead of presenting the frequency depen- 
dence of He it is instructive to plot the dependence of He 
on the applied field magnitude Hq for fixed frequencies. 
Fig 13 shows the Hc{Hq) dependence in the broad range 
(up to 200 Oe) of i^o at the three frequencies, as extracted 
from Fig {To] data. Within this field range transformation 
of minor- into the major-hysteresis loops inevitably takes 
place (Figs. 7|1Q ). Deviation of the experimental points 
in Fig 13 from the present minor loop model (solid lines) 
can naturally be interpreted as a consequence of the lat- 
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FIG. 13. Coercive field He as a function of applied field am- 
plitudes Hq for three measuring frequencies. The plots are 
derived from measurements shown in Fig|10[ Solid lines rep- 
resent the He values calculated from Eq[4] using the same 
T{f,Ho) functional dependence involved in the thick-line fit 
specified in Fig |ll| with the parameters listed in Fig |12| 



ter transformation. At the highest frequency (60 Hz) 
the functional form He = Hc{Hq) is not only linear but 
the simple equality He = Hq holds, up to high values of 
Hq. Such a behavior is fully consistent with the model's 
prediction in the high-frequency {ujt ^ 1) limit (hori- 
zontally aligned loop). For measurements at lower fre- 
quencies, however, Hc{Hq) deviates from linearity and 
tends to saturation already in much smaller values of 
Hq. By comparison with master plots Fig 10 it is clear 



that the Hc{Hq) saturation tendency reflects the stages 
of the minor- to-major loop transformation. Whatever 
the frequency, as long as Hq value grants the whisker- 
free (i.e., the minor loop) regime Hc{Ho) keeps its lin- 
ear slope; for measurements at 60 Hz it stays true up 
to 100 Oe or more (cf. Fig 



however, deviation from the 



13). At lower frequencies. 



inear Hc{Hq) slope sets-in; 



lower the frequency stronger the deviation from linearity. 



Referring to master plots FigjTOJone notes that the quasi- 
saturation level in Hc{Hq) dependence, diminishing with 
frequency, is correlated with the constantly growing size 
of the dynamics- free whisker region in increasing Hq. 



V. DISCUSSION AND CONCLUSION 

Result that the experimental curves become consistent 
with the present model of thermal-activation-enhanced 
Debye relaxations, simply by letting the relaxation time 
r to depend on field and frequency (Fig|l2|, deserves 
some additional comments and considerations. It is par- 
ticularly interesting to discuss the range of applicability 
of this model in the frequency domain by investigating 
the matching of the two coercive fields, dynamic (subject 
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of this work) and the static one. 

Let us first summarize some general knowledge on the 
domain wall dynamics underlying hysteresis loops. Dur- 
ing a hysteresis field cycle the collective spin reversal 
is subject to different dynamical regimes ruled by dif- 
ferent microscopic mechanisms. As a result, within a 
single cycle the domain walls are subject either to wall 
propagation or the domain nucleation^^. Dy namic s of 
these modes involves different relaxation times^^'^ and 
their activation depends not only on the actual mag- 
netic field value (and its history) but also on frequency 
of the applied field (field ramping rate). It has been 
suggested^, for example, that the mechanism of wall 
propagation (with its specific relaxation time) dominate 
over the mechanism of domain nucleation at low frequen- 
cies (and vice versa at high frequency), introducing a 
minimum in the frequency dependence of coercive field 
at the characteristic -cross-over- frequenc}^^^. 

Let us apply this understanding to analysis of our re- 
sults/model. Within a single field cycle of our sinusoidal 
field there is, as discussed above, processes involving very 
different 'instantaneous' relaxation times. These relax- 
ation times, accumulated during a cycle, directly infiu- 
ences properties of our experimental dynamic hysteresis 
(like He) through fundamental relationship for the phase 
shift in our model, (p = arctancjr. We can discuss now 
the observed field/frequency dependence of our effective 
r in more details. The central dissipative process, collec- 
tive spin reversal, can be modeled as the collective spin 
jumps over the energy barrier separating the two neigh- 
boring energy minima. In micromagnetic models^ most 
of magnetic dynamics' features (like the relaxation times) 
evolve from properties of the corresponding energy land- 
scape. The barrier heights, thus the whole landscape, de- 
pend sensitively on the applied magnetic field. The role 
of the applied ac field is therefore twofold: First, it con- 
tinuously redistributes population of the 'up' and 'down' 
domains within the field cycle to keep track with the 
Boltzmann distribution function. Redistribution takes 
place through the dissipative movement of domain walls. 
Second, magnetic field reshapes, through local reduction 
of the barrier heights, the energy landscape inside which 
the collective spin reversal takes place. Field-induced 
reshaping of the energy landscape has a direct conse- 
quence on the effective relaxation times^^ - shallower the 
barriers, shorter the relaxation time. Thus even weak 
fields involved with minor loops reduces the effective r, 
in agreement with our results for rh{Ho). 

In stronger fields, involved with major loops, there are 
different sources of the effective-r reduction. First, the 
domain nucleation processes are now necessarily involved 
contributing with shorter relaxation time in the overall 
phase lag. Second, the frequency and field effects on the 
effective relaxation are actually mixed together: In the 
saturated (whisker-like) part of major loops there is no 
dynamics at all, as the domain walls had been ejected 
out of, now monodomain, sample. Dynamics exists only 
in the non-saturated parts of the loops which, at fixed 



E 0.02 






Hdc(Oe) 

FIG. 14. DC SQUID hysteresis of Co7(Te03)4Br6 at 20 K. 
Field cycle time took approximately 20 hours. 



field frequency, are ramped by the field with the sweeping 
rates proportional to the total field swing. As mentioned 
above, high sweeping rates (equivalent ly, high frequen- 
cies) select out the processes involving fast relaxations. 

As far as the frequency dependence of the effective r is 
concerned it can generally be interpreted along the same 
track- low frequency dynamics relies on processes with 
long relaxation times, typically involved with the pinned 
domain walls^. Important aspect of the frequency de- 
pendent T/(/) is that it naturally bridges our dynamic 
hysteresis loops with 'static' hysteresis. Strictly speak- 
ing, dynamic He vanishes in a true dc limit (cj ^ 0) in 
our model. However, any dc hysteresis can be more re- 
alistically and accurately represented as an ac one but 
taken in the low-frequency limit. For example. Fig. 14 
shows the dc-SQUID hysteresis data at 20 K. In taking 
this hysteresis loop the field cycle took approximately 
20 hours, corresponding to frequency of / = 1.4 x 10~^ 
Hz. The time spent on field cycling (at the order of 
hours) was found infiuencing the shape and the width 
(i.e.. He) of 'dc' hysteresis loop s, as expected in mod- 
els relying on thermal activatiorF*^. However, experi- 
mental value of H^'^{^ 20 Oe, Figjl4|) for mentioned / 
is much bigger than the estimation for Hc{oJ^r) on ba- 
sis of Eq.Q for the corresponding choice of ujt (where 
u = 27r/ and r are calculated as r(/, Hq) = rf{f)rh{Ho) 
with the parameters specified in Fig. 12). Obviously, the 



low- frequency behavior of t/(/) (for / <C 1 Hz) cannot 
be easily extrapolated from studies of dynamic hysteresis 
in the frequency range / > 1 Hz). The reason is a com- 
plicated nonlinear physics of pinned domain walls^^^^^ 
which cannot be, at this stage, easily reduced to a one 
single phenomenological parameter, r. The triangle area 
in FigfTTj indicates a region in the Hc{f) plot character- 
ized by pronounced domain wall pinning effects. In the 
rest of the plot area the thermal-activation-enhanced De- 
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bye relaxation model, presented in this work, obviously 
perfectly interprets the dynamic hysteresis experimental 
data of Co7(Te03)4Br6. 

In conclusion, we studied slowing down of the domain 
wall dynamics below the ferromagnetic transition of a 
new oxo-halide system Co7(Te03)4Br6 by the induction- 
type dynamic hysteresis loops. Convenient quantitative 
matching of the dynamics' time scale with the available 
frequency window enabled detailed studies of minor loops 
in their dependences on frequency and applied field am- 
plitude. The presented model for dynamic hysteresis 
loops relies just on the traditional ingredients of the spin- 
lattice relaxation theory and on the Neel-Brown thermal 



activation. The very low frequency deviations of exper- 
imental results from the model predictions are ascribed 
to specific physics of pinned domain walls, which remains 
outside scope of this work. 
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